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A NOTE ON THE DRAG DUE TO LIFT OF RECTANGULAR 
WINGS OF LOW ASPECT RATIO 
By Edward C. Polhamus 


SUMMARY 


Methods of estimating the induced drag of low-aspect -ratio wings 
are discussed and compared with experiment. The profile drag due to lift 
is also discussed and a method is developed which relates the effect of 
aspect ratio on the profile drag due to lift to an ^effective” two- 
dimensional lift coefficient. A simple expression for this effective 
two-dimensional lift coefficient in tenns of the aspect ratio is derived 
and used to correlate experimental values of profile drag due to lift for 
rectangular wings in the low-aspect -ratio range. A method of utilizing 
two-dimensional test results, by means of this effective two-dimensional 
lift coefficient, to estimate the profile drag due to lift of low-aspect- 
ratio wings is briefly discussed. 


INTRODUCTION 


Considerable work, both experimental and theoretical, has been done 
with regard to the aerodynamic characteristics of low-aspect -ratio wings. 
Comparisons of the experimental and theoretical results have indicated 
that, at least for incompressible flow, the lift and pitching-moment 
characteristics of low-aspect -ratio wings can be estimated with reason- 
able accuracy. (See, for example, ref. 1.) The estimation of the drag 
due to lift, however, has been hampered by what appears to be a large 
effect of aspect ratio on the variation of profile drag with lift coef- 
ficient. The purpose of the present paper, therefore, is to attempt to 
determine the effect of aspect ratio on the variation of the profile drag 
with lift coefficient . 


SYMBOLS 




section lift coefficient. 


Section lift 
qc 


^effective" section lift coefficient, 


"Effective” section lift 


qc 
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ACd 

ACOp 

CDi 

Cn 


Cs 
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q 

p 

V 


C 
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total lift coefficient^ 


Lift 

qS 


portion of total lift coefficient associated with induced 
Induced-camber lift 


camber^ 


qS 


section drag coefficient^ 


Section drag 
qc 


total drag coefficient^ 


Drag 


drag coefficient at zero lift 


drag -due -to -lift coefficient^ ■ *^Do 


coefficient of profile drag due to lift, ACp. - C 


D ‘ ^D^ 


induced-drag coefficient. Induced dra.g 

qS 


normal-force coefficient. 


Normal force 
qS 


re suit ant -force coefficient, Resultai ^ force 

qS 


section suction-force coefficient. 


Section suction force 
qc 


suction-force coefficient, S uct J^ pn force 

qS 


lift -curve slope 


dynamic pressure, pV^/2, Ib/sq ft 
mass density of air, slugs/cu ft 
free-stream velocity, ft /sec 
wing chord, ft 
wing area, sq ft 
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A 

E 

b 

a 

X 

AP 

P 

Po 

P 


aspect ratio, b^/S 

Jones' edge-velocity correction factor 
wing span, ft 

geometric angle of attack, radians 
longitudinal distance from wing leading edge , ft 

pres sure -difference coefficient, Pupper " Plower 

local static -pres sure coefficient, ^ 

free-stream static pressure, Ib/sq ft 
local static pressure, Ib/sq ft 

ANALYSIS AND DISCUSSION 
Theoretical Induced Drag 


The lifting-line theory for induced drag developed by Prandtl 
(ref. 2) states that for an elliptical spanwise load distribution the 
induced drag is given by 



( 1 ) 


The lifting-line theory, of course, neglects any effect of the vortex 
sheet on the chordwise loading due to induced curvature of the stream- 
lines, and therefore it is not applicable to low-aspect -ratio wings with 
regard to lift and moment. Lifting-surface solutions must be used to 
account for this effect of the vortex sheet. With regard to the induced 
drag, however, Munk has pointed out (ref. 3) that the total induced drag 
is independent of any staggering of the lift in the stream direction and, 
therefore, the Prandtl lifting-line solution, with regard to induced drag, 
is applicable to all types of wings, provided the spanwise loading is 
elliptical. The theory assiimes, however, that the vortex sheet lies in 
the plane of the wing and therefore it is strictly correct only for van- 
ishingly small angles of attack. In an attempt to account for the effect 
of angle of attack. Mangier, following a proposal by A. Betz, assumes 
that the displacement of the vortex sheet can be treated in a manner 
similar to that used to predict the effect of end plates (see ref. 4 ). 
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Tnc xieight of the imaginary end plates increases with angle of attack, 
and thus the effective aspect ratio is increased. The resulting 
expression for the induced drag is 


CDi = 


Cl" 


jt(A + a) 


( 2 ) 


with a in radians. (it is interesting to note that an essentially 
equivalent result has been obtained by Weinig (ref. 5) with a momentum 
approach.) Inasmuch as the difference between equations (l) and (2) is 
appreciable only for aspect ratios of about 1.0 or less, an aspect ratio 
of 0.25 has been chosen for a comparison of the two theories with experi- 
mental results. This comparison is shown in figure 1 for a rectangular 
wing, the tests of which are described in reference 6. Inasmuch as the 
leading-edge radius of this wing is relatively large (approximately 0.021c), 
potential flow would be expected up to moderate lift coefficients and would 
afford a good check on the potential-flow theories. It should be pointed 
out that in the interest of accuracy in the low-lift range, the data pre- 
sented in figure 1 were obtained from the original tabulated data and not 
from the plot presented in reference 6. The results indicate that the 
drag due to lift of this wing is considerably lower than the theoretical 
induced drag as given by Cj^/Tik and is in good agreement with that 

given by Cj^ V«(A + a) for lift coefficients up to about 0.20. At the 
higher lift coefficients, apparently nonpotential flow effects cause the 
experimental drag due to lift to rise more rapidly than the theoretical 
curve for induced drag (eq. (2)). This difference in drag is considered, 
in the following sections, as profile drag due to lift. 


Effect of Aspect Ratio on Profile Drag Due to Lift 

In the preceding section, it has been shown that in the low- lift 
range, where the profile drag is relatively independent of lift coeffi- 
cient, the drag due to lift can be predicted fairly accurately by equa- 
tion (2). The problem remaining in the prediction of drag due to lift, 
therefore, is that of the variation of the profile drag with lift. In 
figure 2, the increment of profile drag due to lift determined as 

the difference between the measured drag due to lift and the theoretical 
induced drag (given by eq. (2)), is presented as a function of lift coef- 
ficient for several aspect ratios. The data were obtained from refer- 
ence k and the wings were of rectangular plan form (A = 0.5^ 1^ and 2) 
with NACA 0012 root sections and NACA 0009 tip sections. The wings were 
tested at a Reynolds number of about 106. The results indicate that the 
profile drag due to lift is dependent to a large extent upon the aspect 
ratio of the wing, with extremely high values occurring for the low-aspect 
ratio wing. This, of course, indicates that the profile drag for a 


NACA TN 5524 


5 


three-dimensional ving at a given lift coefficient cannot he determined 
from corresponding two-dimensional tests at that lift coefficient. In 
figure 5 ^ a similar comparison is made for wings having thinner airfoil 

sections (NACA 65 AOO 6 ) at a Reynolds number of approximately 2 X 10^. 

These data were obtained from semi span -model tests (unpublished) conducted 
at the Ames Aeronautical Laboratory of the NACA. Although the profile 
drag due to lift ACpp is considerably greater for these wings than for 

the thick wings of figure 2, the effect of aspect ratio is similar. The 
reason for the higher drag on these thin wings is that the adverse pres- 
sure gradient at the leading edge increases with decreasing leading-edge 
radius and thereby causes more severe separation on the thin wings. This 
phenomenon also accounts for the fact that the data of figure 1 are not 
consistent with those of either figure 2 or figure 5 * 

The large effect of aspect ratio on the profile drag may be associated 
with the fact that, as the aspect ratio decreases, the chordwise pressure 
gradient due to angle of attack in the vicinity of the leading edge 
increases as a result of the induced camber (streamline curvature) associ- 
ated with the chordwise variation of the induced downwash. (See, for 
example, ref. 7*) This increase in the adverse pressure gradient, through 
its effect on separation at the leading edge, would be expected to cause 
an increase in the profile drag. This effect of aspect ratio on the 
chordwise pressure gradient is illustrated in figure 4. In the upper 
right-hand corner of this figure the theoretical two-dimensional chord- 
wise pressure distribution due to angle of attack over an NACA 0012 air- 
foil (ref. 8 ) is compared with the theoretical chordwise loading at the 
root section of a rectangular wing of aspect ratio 0.5 that also has an 
NACA 0012 airfoil section.^ The change in the chordwise pressure distri- 
bution with decreasing aspect ratio is clearly evident in figure 4 and 
results, as pointed out previously, in larger pressure gradients in the 
vicinity of the leading edge. These differences in chordwise pressure 
distributions are a result of the induced camber (streamline curvature) 
associated with the chordwise variation of induced downwash for the wing 
with aspect ratio of 0 . 5 . This induced camber corresponds to a negative 
geometric camber and increases in magnitude with decreasing aspect ratio. 

In order to maintain a given lift coefficient the loss due to the induced 

^Available three-dimensional calculations are for wings of zero 
thickness and result in infinite pressures at the leading edge. For this 
comparison, however, it is more satisfactory to utilize wings that have 
finite pressures at the nose and, in addition, have the same airfoil sec- 
tion as was used for the experimental results presented in figure 2 . 

Since from reference 9 can be shown that the chordwise distribution 
of induced camber at the root of a rectangular wing is similar to the 
camber distribution of the NACA 65 mean line, the corresponding chordwise 
pressures (p. 9 I of ref. 8 ) corrected to the proper camber load with the 
aid of figure 6 were combined with the angle -of -attack loading of an 
NACA 0012 airfoil (p. 7 I of ref. 8) to obtain the results shown for A = 0 . 5 . 
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camber, which is a maximum at the midchord and decreases to zero at the 
leading and trailing edges, must be compensated for by an increase in 
angle of attack. The angle-of -attack load, 'however, is a maximum at the 
leading edge and therefore, although the loss of lift is compensateci for, 
the distribution along the chord is altered as illustrated in figure k. 

In the lower part of this figure the chordwise pressure distribution near 
the nose has been "blown up" to facilitate a comparison. It can be seen 
that the magnitudes of the pressure near the leading edge of the wing 
with aspect ratio of 0.5 (solid curve) are considerably greater than 
those for the two-dimensional wing (long-dash curve) at the same lift 
coefficient, and that the lift coefficient of the two-dimensional wing 
would have to be approximately twice that of the wing with aspect ratio 
of 0.5 (long- and- short -dash curve) to develop the same pressures near the 
nose. This increase in lift coefficient is equal to the negative lift 
associated with the induced camber. It would therefore appear that 
boundary-layer and separation characteristics of the wing with aspect 
ratio of 0.5 inay more nearly correspond to those of the two-dimensional 
airfoil operating at twice the lift coefficient. This may account for 
the large effect of aspect ratio on the variation of profile drag with 
lift coefficient shown in figures 2 and 3 - 

Before attempting to correlate the effect of aspect ratio on the 
profile-drag variation with lift on the basis of the above concept, it 
is interesting to note that this concept accounts for the fact, mentioned 
previously, that lifting-line and lifting-surface theories result in the 
same induced drag despite the fact that lifting-surface solutions result 
in a more rearward inclination of the normal force (lower Cl^^) . In 
order to make up for this increased reai’ward inclination, the suction 
force at the nose is required to be greater by lifting- surface concepts 
than by lifting-line concepts. The values of the suction-force param- 

eter C 3 /CL required to incline the resultant force at the correct 

angle relative to the free stream were determined from the lifting-line 
and lifting-surface theories (see appendix) and are compared in figure 
It will be noted that at low aspect ratios considerably greater suction 
is required than is given by lifting-line or two-dimensional theory but 
that when the suction is calculated by a two-dimensional theory (eq. (A9) ) 
that incorporates a lift coefficient increased by ACj^, which corresponds 

to the loss due to induced camber, good agreement is obtained. The reader 
is referred to the appendix for a more detailed discussion of figure 5 * 

It now appears possible that the increment of profile drag due to 
lift may be independent of aspect ratio for a given value of the effec- 
tive two-dimensional lift coefficient 
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where AC^ is equal to the amount of induced-camber load. In the appen- 
dix^ several methods of determining the induced-camber load are discussed 
and it is shown that to a very good approximation (for unswept wings) 



AC 


A plot of 1 + — is presented in figure 6 and it will 

low aspect ratios the induced -camber load becomes rather 
son with the total load. 


be noted that at 
large in compari- 


Figures 7 s-nd 8 have been prepared in order to evaluate the degree 
to which the use of this effective two-dimensional lift coefficient will 
eliminate the effect of aspect ratio on the profile drag due to lift. In 
figure 7^ data of figure 2 for profile drag due to lift are replotted 

as a function of the effective two-dimensional lift coefficient 



The results indicate (fig. 7) that the effect of aspect ratio shown in 
figure 2 can be eliminated to a large extent by use of this effective 
two-dimensional lift coefficient. In figure 8, the same comparison is 
made by use of the data for the thin wings presented in figure 5^ and^ 
although the values of ACp^ are considerably higher than those for the 

thick wings (fig. 7), the effect of aspect ratio again is largely elimi- 
nated. It is probable that the degree of correlation indicated is as 
good as could be expected, inasmuch as effects such as those caused by 
possible changes in span loading with angle of attack, tip profile shape, 
and so forth, although relatively small in comparison with the profile 
drag due to lift, have not been accounted for in the determination of the 
induced drag. 

With regard to the application of two-dimensional data in the estima- 
tion of three-dimensional drag due to lift, the preceding results indi- 
cate that two-dimensional results for a lift coefficient that is higher 

^by the factor -1^ than the lift coefficient of the three- 

dimensional wing in question must be used. 
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CONCLUDING REMARKS 


An important problem in connection vlth the prediction of the drag 
due to lift of lov-aspect -ratio wings involves the determination of the 
effect of aspect ratio on the variation of profile drag with lift coef- 
ficient. In this connection the possibility of correlating the values 
of profile drag due to lift for wings of various aspect ratios by the 
use of an "effective** two-dimensional lift coefficient has been shown 
both theoretically and experimentally. It is shown that, to a good 
approximation, this effective two-dimensional lift coefficient can be 


and A is the aspect ratio. When the profile drag for several aspect 

ratios was plotted as a function of ^ effect of aspect 

ratio was eliminated to a large extent. This, of course, implies that 
when utilizing two-dimensional data to determine the three-dimensional 
profile drag, the two-dimensional profile drag corresponding to a higher 
lift coefficient must be used. Other implications are that effects that 
occur in two-dimensional flow (such as Reynolds number, profile shape, 
and so forth) will occur at progressively lower lift coefficients and be 
more severe as the aspect ratio is reduced. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 28, 195^* 


expressed as 



where C-^ is the total lift coefficient 
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APPENDIX 


DERIVATION OF AN "EFFECTIVE" TWO-DIMENSIONAL 


LIFT COEFFICIENT 


As pointed out in the body of the present paper, lifting-line and 
lifting-surface solutions result in the same expression for the induced 


theory is considerably lower, especially at low aspect ratios. This 
means that, despite the fact that the inclination of the normal force 
relative to the free stream given by lifting- surface theory is greater 
than that given by lifting-line theory, the inclination of the resultant 
force is the same by both theories. In order to produce the same incli- 
nation, the suction force at the leading edge, which rotates the resultant 
ahead of the normal, must be greater than that given by lifting-line 
theory. This is illustrated in figure 4 where the leading-edge suction 

parameter Cg/Cj^ as determined by the various theories is compared. 

For a two-dimensional wing^ since the suction force must, when combined 
with the normal force, produce a resultant lift force normal to the free 
stream: 


drag 



although the lift-curve slope given by lifting- surface 



(Al) 


Therefore 


cj2 2rt 


(A2) 


For a three-dimensional wing the suction must be such as to incline 
the resultant force rearward by the angle Therefore 
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By lifting-line theory. 


Cl(A + 2) 
2rtA 


Substituting the above equation into equation (A5) yields 


Cs = Cl 


E 


Cl(A + 2) 


2jtA 



which reduces to 



This equation is independent of aspect ratio and identical to the two- 
dimensional value which is given by equation (A2) , and is shown in fig- 
ure 4 as the long-and-short-dash curve. However ^ when the angle of 
attack a determined by a lifting-surface solution (ref. 10, for example) 
is substituted in equation (A5), it is found that the suction is greater 
than indicated by lifting-line theory and dependent upon the aspect ratio 
as shown by the solid line in figure 4. In order to extend the results 
to aspect ratios approaching zero, the slender-wing solution of Jones 


(ref. 11), a 
in 


2C 

may be substituted into equation (A5) . This results 


P 2 kA 


(A4) 


and the resulting cui^e is presented in figure 4 as a long-dash line. 
Although the theory of reference 11 is for vanishingly small aspect 
ratios, in general it gives accurate results for the lift-curve slope 
up to aspect ratios of about 0.5* The dashed line in figure 5 is there- 
fore extended up to A = 0.5 and it appears to be a reasonable extension 
of the results obtained with the aid of reference 10. 

It would be interesting to see now if the increased pressures at 
the leading edge due to the induced camber mentioned in the body of the 
paper will account for the difference in suction required by lifting-line 
and lift ing- surface theory. Now, equation (A1) gives the suction as 
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This equation is correct, however, only when the induced angle is zero 
as in two-dimensional flow or is constant along the chord as is assumed 
by three-dimensional lifting-line theory. For the actual three-dimensional 
case, however, where the induced angle is not constant along the chord, 
the pressures at the nose correspond to those of a two-dimensional wing 
operating at a higher lift coefficient as illustrated in figure 4. This 
effective two-dimensional lift coefficient can be expressed as 


where ACj^ is the magnitude of the induced-camber load. The problem 
now reduces to the determination of the induced-camber load. 

In reference 12 Jones corrected the Prandtl lifting-line equation 
for the lift-curve slope of elliptical wings for the effect of the reduc- 
tion in velocity around the aft semiperimeter of the wing. This reduc- 
tion, which is due to the fact that the semiperimeter is greater than 
the wing span, results in less circulation being required to satisfy the 
Kutta condition and thereby less lift. This "edge velocity" correction 
results in the following equation (for a section lift-curve slope of 2it) 



and equation (Al) for the suction can be rewritten as 



or 




(A5) 



(a6) 


where E is the ratio of the semiperimeter to the span of the elliptical 
wing and is given by the following expression 
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Now it would appear that the only difference between the results obtained 
by equation (a 6) and the exact lifting- surface solutions is the induced 
camber associated with the variation of downwash along the chord. 
Therefore 



(%)e 


(A8) 


where C]^ is the lifting-surface solution. 

Substitution of equation (A8) into equation (A5) results in 



(A9) 


Now by utilizing the lifting-surface results for elliptical wings 
(ref. 10) for Ci^ and equation (a 6) and (AT) for (Cr ') , the suction 

force has "been calculated hy equation (A9) and the results com- 

pared with the exact results in figure 5- The agreement obtained appears 
to substantiate the assumption that the flow about the leading edge of a 
finite-span wing is equivalent to that about a two-dimensional wing at an 
effective lift coefficient that is greater than the actual finite-span 
lift coefficient by an amount which can be determined by use of equa- 
tion (A8) . 


An alternate method of determining the effective two-dimensional 


lift coefficient 



which does not require that the induced 


camber be known is as follows. If it is assumed that the flow is non- 
viscous and potential^ except that no suction force is developed at the 
leading edge^ so that the resultant force is normal to the chords the 
total drag for both two- and three-dimensional wings is equal to the lift 
times the tangent of the angle of attack: 


C|^ = c^ tan a 



2rr 


(AlO) 
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and 


Cjj = tan a 



(All) 


However, the three-dimensional drag is also equal to the corre- 
sponding two-dimensional drag plus the induced drag 


Cd = Cd -H 




Therefore, from equation (All), 



Solving for c^ gives 



(A12) 


(A13) 


(Aik) 


Now if the lifting- line equation for Cj. , 




2nA 
A + 2 


is substituted into equation (A14) it can be shown that c^ = How- 

ever, if lifting-surface solutions are substituted for Ci^ in equa- 
tion (A14), it can be shown that in order for the two-dimensional drag 


I 


Cl 


plus the induced drag 4- to equal the three-dimensional drag 


2jt 


jtA 


Cl /Ci^ section-lift coefficient must be greater than the three- 


dimensional lift coefficient C]^. This is the same conclusion reached 
previously when the induced camber was considered; a comparison of the 
magnitudes of the increases given by the two methods is made in the fol- 
lowing paragraphs . 
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From equation (A14), 


1 



(A15) 


where ACj^ is the increase over the three-dimensional value of C^. 

From equation (a 8), which was determined from a consideration of the 
induced camber. 


1 


+ 


Cl 


1 + 2it 


1 

Cla 




(ai6) 


The results of these two methods are compared in the following table 


where Cj;^ is determined from reference 10 and 
tions (a 6) and (A7). 


(^^)e 


from equa- 


A 

% 


ACt 

from — 



eq. (A15) 

eq. (A16) 

6.57 

4.55 

4.6l 

1.05 

1.02 

2.55 

2.99 

3.15 

1.15 

1.11 

1.272 

1.82 

2.00 

1.57 

1.51 

.637 

• 99 

1.15 

1.79 

1.79 


It will be noted that fairly good agreement is obtained between the 
two methods and, inasmuch as equation (A15) is dependent upon results from 
only one theory whereas equation (A16) is dependent on the results of two 
theories, it is felt that equation (A15) niay be somewhat more reliable 
and it is, of course, more convenient to use. Now if the simple but very 
accurate expression for developed by Helmbold (see refs. 15 and 1^) , 

Ct = — + 2 , is substituted into equation (A15) , the following is 

obtained: 
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(A17) 


Values obtained from this expression are compared in figure 6 with 
the results obtained from equation (A15) by using Cj_^ from reference 10 
and good agreement is indicated. 
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Figure 1.- Comparison of experimental drag due to lift with two 

theories for induced drag. Symbols represent experimental data from 
reference 6. 



^Cop = A Co - 
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Figure 2.- Effect of wing aspect ratio on the profile drag due to lift 

NACA 0012 root section and NACA 0009 tip section. Data are from refer- 
ence 4. 
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Flgxire 3 .- Effect of wing aspect ratio on the profile drag due to lift 

NACA 65 AOO 6 airfoil section. Data are from tests conducted at the Ames 
Aeronautical Laboratory. 
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Fig-ure 



Calculated effect of aspect ratio on the chordwlse load distri- 
bution due to angle of attack. NACA 0012 airfoil section. 
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Figure 5»“ Variation of leading-edge suction parameter with aspect ratio 

as predicted by several theories. 
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Figure 6.- Effect of aspect ratio on the induced-camber parameter 1 + — 


as predicted by two methods 



Figure 7.- Correlation of the profile drag due to lift of figure 2 with an 

effective" two-dimensional lift-coefficient parameter, ro 
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Figure 8.- Correlation of the profile drag due to lift of figure 5 with 
an "effective" two-dimensional lift-coefficient parameter. 
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